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Abstract 

Using a recursive algorithm to solve the renormalization group equations 
of = 1 QCD (DGLAP), we describe the most general supersymmetric evo- 
lution of the parton distributions. The analysis involves the regular DGLAP 
evolution, a partial supersymmetric intermediate evolution and a final super- 
symmetric evolution combined at various matching scales. We use a model 
in which supersymmetric distributions are radiatively generated at each susy 
threshold due to the mixing of the QCD anomalous dimensions with the = 1 
sector. Various types of matching conditions are considered, reflecting partially 
broken and exact supersymmetric scenarios with a wide range of susy-breaking 
mass parameters. Numerical results for all the distributions are presented. 
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1 Introduction 



The search for viable supersymmetric extensions of the Standard Model has acquired 
a new momentum with the coming into operation of the Large Hadron Collider at 
Cern, scheduled for the near future. A great deal of work is currently undertaken 
in trying to map the various scenarios that supersymmetric models predicts in this 
unexplored energy range, by providing estimates for the various channels which may 
become available at this new scale. Being the symmetry broken at such energy, the 
several parameters describing the breaking make its description a complex task. In 
this respect, there are different avenues that these studies can take. For instance, 
one possibility is to investigate supersymmetric corrections affecting the initial state, 
prior to the hard scattering phase; or to analize the opening of heavy supersymmet- 
ric channels in the intermediate/final state of the collisions; or, finally, to investigate 
a combination of both. These avenues are all - although not equally - justified by 
current experimental lower bounds on the mass of the supersymmetric partners. Su- 
persymmetric scaling violations, which are those considered in this work, have to do 
with the first type of searches. They have also to do with a regime of the theory 
in which the symmetry is - to some extent - restored, and all the channels become, 
effectively, massless. We will elaborate on them in some detail in the following sec- 
tions. 

We have presented in a previous paper results for the evolution of regular and 
supersymmetric parton distributions within a scenario characterized by broken su- 
persymmetry and a decoupled squark. In this work we intend to examine this subject 
from a different perspective. We analize an exact supersymmetric phase of the the- 
ory and study the structure of the distribution of quarks and gluons and of their 
supersymmetric partners, gluinos and squarks, under these conditions. This scenario 
is, obviously, less realistic than a scenario characterized by a broken susy, but not 
less interesting neverthless. In particular, we will present some comparisons between 
the regular QCD scenario, various types of susy breaking scenarios and the exact 
supersymmetric scenario. 

The aim of this analysis is to quantify these effects with a reasonable accuracy, in 
such a way that the results can be used as a guideline for experimental searches in the 
future. We build on previous old work of Kounnas and Ross who calculated the 
leading order anomalous dimensions for the evolution of supersymmetric QCD (Susy 
DGLAP) and analized the first 2 moments of the distributions. Our phenomenological 
analysis is updated to current estimates of the parameters of the parton distributions. 
It includes all the moments, since it is an exact numerical solution of the evolution 
performed by standard iteration of the convolution integrals, recast in the form of 
a recursion relation. In the section that follows, we will use the notation AP to 
denote the standard DGLAP evolution. The acronym SDGLAP or, briefly, SAP, 
will denote the partial supersymmetric DGLAP evolution with coupled gluinos and 
decoupled squarks. Finally, the acronym ESDGLAP or ESAP will denote the exact 
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supersymmetric DGLAP for iV = 1 QCD with coupled gluinos and squarks. 



2 The Evolution Equations of SQCD with Exact 
Super symmetry (ESAP) 

We refer to previous work of us for a detailed description of the algorithm that we 
are going to use in this analysis ^j. We introduce some definitions just in order to 
make our discussion self contained. 

Similarly to ordinary QCD, we define singlet and non-singlet distributions 
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The evolution equations can be separated in two non-singlet sectors and a singlet 
one. The non-singlet are 
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where the non singlet (NS) kernel is given by 
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and the singlet, which mix qy and qy with the gluons and the gluinos are 
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There are simple ways to calculate the kernel of the supersymmetric evolution 
by a simple extension of the usual methods. The changes are primarily due to color 
factors. There are also some basic supersymmetric relations which have to be satisfied 
that will be analized below. They are generally broken in the case of decoupling. We 
recall that the supersymmetric versions of the (3 functions are given at 1-loop level 
by 
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with ur = riL = Tip = nf/2. The running of the coupling is given by 
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We use the ansatz discussed in P| which amounts to take standard iterates of the 
convolution products in a large number to solve the equations, implemented in the 
form of recursion relations fixed at runtime. We start applying the method to the 
non-singlet case and then move to the singlet. 



3 Models of evolution 

We are going to discuss simple models of the supersymmetric evolution which incorpo- 
rate several phases: the N=0 or regular QCD phase, with fully decoupled superpart- 
ners, described by the usual DGLAP equations; the N=l phase, characterized by a 
decoupling of some of the superpartners (the squarks) from the remaining evolution, 
and, finally, the exact supersymmetric phase in which all the fields are effectively 
massless and are evolved simultaneously. The last two phases can be obtained as 
intermediate and final stages of an evolution with a gluino mixing. 

Beside the initial scale Qi, at which we start the evolution, the beginning and 
the end of the intermediate phase (SAP) will be denoted by two additional scales, 
denoted by Qs and Qes- At some point, we will choose to vary the location of these 
two scales, and analize the impact of the evolution according to all possible scenarios 
which are implied by these choices. 

Some of these scenarios will be unrealistic, while others are more supported by 
the current experimental lower bounds on the masses of the superpartners. All the 
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cases that will be discussed below have the objective to illustrate in a sufficiently 
detailed way the main features of the susy evolution. 

We should mention, if not obvious, that the intermediate and final phases in 
the evolution (evolved by the SAP and the ESAP respectively) are all affected by 
"asymmetric" boundary conditions, which are characterized by some densities of the 
superpartners set to zero at a given intermediate threshold. 

This reasoning, as explained in our previous work, is in line with current ap- 
proaches to the analysis of the QCD evolution in which the parton densities can be 
generated by radiative collinear emissions from low scale distributions. We recall that 
the role of scaling violations and of the renormalization group all-together is, in this 
context, to simply dress light-cone matrix elements by logarithmic enhancements. 

It is convenient to introduce a general notation for all the kernels that will be 
considered. We embed both the singlet AP and SAP kernels into 4x4 matrices 
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while the non singlet DGLAP kernel is rewritten in a 2 x 2 form 



pNS 
^AP 



Pm 




(10) 



For convenience we introduce the following notations. A region in the RG 
evolution - regular, partially decoupled or exact supersymmetric - is described by 
an array {Qi,Qf), where Qi is the initial evolution scale and Qf the final evolu- 
tion scale. The type of evolution (AP,SAP,ESAP) is indicated by a correspond- 
ing suffix. For instance (2, 100)ap denotes a regular DGLAP evolution with initial 
scale Qi = 2 GeV and a final scale of Q/ = lOO GeV. Similarly, an evolution of 
the form (2, 100)ap - (100,400)5ap - (400, 1000)s5Ap describes a matching of aU 
the three evolutions at the intermediate scales Qf,AP = Qi,SAP = 100 GeV and 
Qf,SAP = Qi,ESAP = 400 GeV. The masses of the intermediate scales are denoted 
by an array with 2 entries {m2\,m2q)- In the example presented above, these two 
scales are the matching scales for the SAP and ESAP evolutions Qi^sAP ~ '>Ti2X and 
Qi,ESAP = triq. As usual, we adopt a step approximation in the running, accord- 
ing to which we step into a new region right after crossing the corresponding mass 
threshold. We recall that we use a general logarithmic ansatz for the structure of the 
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logarithmic contributions in which the momentum dependence of the scahng viola- 
tions is parametrized by the actual running coupling. In we have elaborated on 
this point in some detail. We set the expansion 
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and invoke an equality between logarithmic powers at the left-hand-side and at the 
right-hand-side of the equations to get recursion relations for the coefficients A-I'nix)- 
TO be specific, let's consider the most general sequence of evolutions (AP-SAP- 
ESAP) described by the arrays {Qi,Qf)AP - {Qi - Qf)sAP - {Qi,Qf)ESAP, with 
Qf,AP = Qi,SAP = and Qf^sAP = Qi,ESAP = ^q- In this (general) case the 
solution is built by sewing the three regions as 







JQI 

+ rd\ogQ'P^ip{x,a'{Q^))® 

+ r'rf log g2p|'/^p(x,a^^(g2))® 

J mi. 



gy(x,g2) 


qv{x,Q^' 



q^{x,Q' 



in the non singlet and 







(12) 



r G{x,Q}) 1 




r G{x,Q}) 1 


A(x,g2) 







g+(x,g2) 




q+ix,Q}) 


. r(x,g2) . 








+ 



2 

r^rflog Q'P^^{x,a{Q')) 
Jo?. 



+ 



d\ogQ''P^ip{x,a'{Q''))(S> 



+ 



log g2p|'/^p(x,a^^(g2))® 



G{x,Q}) 
A(x,g2) 
g+(x,g2) 


G(x,g2) 
A(x,g2) 

g+(x,g2) 
g+(x,g2) 



G{x,Q}) 


g+(x,g2) 




(13) 



for the singlet solution. The zero entries in the arrays for some of the distributions are 
due to the boundary conditions, since all the supersymmetric partners are generated, 
in this model, by the evolution. The general structure of the algorithms that solves 
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these equations is summarized below. We start from the non singlet. We define 

q^^{x,Q'^) = {qv{x,Q'^),qv{x,Q'^)f and An{x) = (A^v',^!^^ and introduce the 
ansatz 
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where no is an integer at which we stop the iteration. Usually ranges between 30 and 
40. The first coefficient of the recursion is determined by the initial condition 
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The recursion relations are given by 
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The solution in the first (DGLAP) region at the first macthing scale m2x is given by 
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At the second stage the (partial) supersymmetric coefficients are given by (S is a 
short form of SAP) 
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We construct the boundary condition for the next stage of the evolution using the 
intermediate solution 

evaluated at the next threshold rriq. 
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The final solution is constructed using the recursion relations 

A^^{x) = U{x) (g)q{x,mq) 



(22) 



The final solution is written as 
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4 Solving the non singlet equations 

In the case of exact supersymmetry, the arguments and the strategy presented in ||^ 
simplify, since we neglect all the intermediate scales and the evolution -starting from 
the lowest scale- is assumed to be supersymmetric. The boundary contions at the 
start of the evolution, however, are not. We implement the algorithm as it has been 
formulated in our previous work, with due modifications given the different structure 
of the evolution kernels. 

The recursion relations needed in the implementation are defined in terms of 
recursive coefficients A-^„(x) with n = 0, 1, 2, ... 

This equation can be made explicit by separating the more singulare terms of the 
recursion relations from the rest 
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5 The Singlet Equation 



The singlet equations are treated in a similar way. The edge-point contributions {x = 
1) appear to be different from the usual QCD expressions for those splitting functions 
that are part of the regular QCD evolution Pgg and Pgg due to the supersymmetry 
relations 

Pgg + Pxg = PgX + PxX (26) 

and by 

Pqg ~l~ Pxq Pgs ~l~ Pxs 
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6 AP-ESAP evolution 

As a starting point, we discuss a model of the evolution in which supersymmetry is 
switched-on right above the ordinary QCD evolution, without introducing an inter- 
mediate region of partial decoupling. This is equivalent to choose mx = rrtq. This 
example may serve as an illustration of the impact of a full (or exact) supersymmetric 
evolution on top of the regular QCD evolution. Initial conditions for the ESAP evo- 
lution are characterized by vanishing densities for all the superpartners at the scale 
Qes where the ESAP evolution starts, and by (leading order) evolved AP distribu- 
tions of ordinary QCD. We plot only distributions summed over all the flavours of 
quarks and squarks such as xg+(a;), xq^{x) and we fix the number of flavours n/ = 4. 

The valence quark distributions qv{x, Qg) and gluon distributions G(x, Ql) at the 
input scale Qq, taken from the CTEQ3M parametrization 

Aox^^{l-x)^^{l + A^x^^). (32) 
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/ N ,u(x) + d(x) , , , 

xs{x) = 0.5x( ^ ^ (33) 

and a vanishing anti-strange contribution at the input. Fig.l shows the shapes of the 
initial CTEQ distributions at an initial scale of = 2 GeV. 



7 AP-ESAP Evolution 

We show in Fig. |1] the shape of the distributions at the lower scale of 3 GeV. They 
are the non-singlet xuv{x), the singlet combination xq^{x) and the gluon xG{x). We 
will be using these initial shapes all along in our analysis. In Figs. and ^ we 
show comparisons between the AP evolution and the combined AP-ESAP evolution. 
As we have already mentioned, this model of the evolution is realistic if all the susy 
partners are close in mass. We start focusing our attention on the regular QCD 
distributions, now evolved in an mixed (regular and supersymmetric) setting. We 
take the matching parameter m2\ to be of 20 GeV, corresponding to a light gluino 
(and to a light scalar quark). The final evolution scale is Q/ = 500 GeV. In Fig. |^ 
we show the evolution of the gluon density in the AP case versus the AP-ESAP case. 
The initial and final scales of the AP evolution are taken to be the same (Qo = 2 
GeV and Qf =500 GeV). Scaling violations are sizeable both for gluons, for valence 
quarks (see Fig.^) and for singlet quarks. As shown in Fig. ^ the singlet squark 
density becomes significant at smaller x. In this figure we also compare in size this 
distribution to the singlet quark one. This distribution is down approximately by 
a factor of 10 compared to that of regular quarks. The decrease at larger-x of this 
supersymmetric distribution is also much faster. 

A similar pattern is shown in Fig. ^ in this case for the gluino density. The 
figure shows the dependence of the density on the matching scale rriq. The gluino 
distributions grows sharply at small x. As expected, as we increase the mass of 
the scalar quark the gluino distribution flattens. In Fig. ^ we compare the valence 
quark distribution for two values of the scalar quark mass in the evolution {rriq = 20 
GeV and 200 GeV respectively). The final evolution scale is fixed to 500 GeV. 
In the same plot we show for comparison the non singlet scalar quark distribution 
(for a mass rriq = 200 GeV). The variation in shape of this distribution due to the 
presence of a supersymmetric threshold in the evolution is comparable to the typical 
scaling violations of ordinary QCD induced by a change of the factorization scale. 
In Fig. ^we plot the non-singlet squark distribution (denoted as sqns) for a varying 
squark mass. The distribution gets lowered drastically as the mass of the scalar quark 
increases. The shape of this distribution, which is generated radiatively -starting from 
an initial scale rriq- is similar to the usual non singlet quark distribution, but rescaled 
by a factor approximately estimated to be of 5/100. Fig. ^ shows the singlet quark 
density for 3 values of m^. The distributions are more pronounced at smaller x values 
and become lowered as the scalar quark mass increases. This is expected, since, for a 
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given supersymmetric evolution interval {Qint, Qf)ESAP, the supersymmetric interval 
gets smaller as Qint —^Qf and therefore, the logarithmic enhancements are reduced. 

We also study the dependence of the evolution on the final evolution scale, keeping 
fixed the matching scale iriq. Figs. ^ and 10 show the dependence of xG{x) from the 
final evolution scale int the AP-ESAP evolution. We have chosen final evolution 
scales of 1,2, and 5 TeV respectively for both plots. The matching scales in the two 
plots are different. We have set iriq = 20 GeV and 200 GeV respectively. The plots 
show a similar pattern and are hardly distinguishable. 

In Fig. 11 we illustrate the dependence of the gluino density on the final evolution 
scales, chosen to be of 1, 2 and 5 TeV. As we increase the final scale Qf for a fixed 
rriq, the gluino distribution becomes steeper. Figs. and |l^ have been included to 



show the dependence of the the non-singlet distributions (valence) for squarks and 
quarks on the final evolution scale Qf. The dependence is slightly more pronounced 
for squarks than for valence quarks. Finally, Fig. |l^ shows the shapes of the gluino, 
gluon and non-singlet squark distributions when the final scale is very large. We have 
taken Qf = 10^ GeV, corresponding to a very energetic collision. In this figure all 
the distributions -except one- have been obtained by an AP-ESAP evolution. One 
of the distributions is evolved using a general supersymmetric run (AP-SAP-ESAP) 
with an intermediate mass gluino (20 GeV) and a slightly heavier squark. 



8 The General SUSY evolution 

The most general supersymmetric evolution is obtained by matching three regions: 
the AP region, the SAP region and the ESAP region. This is also a realistic ap- 
proximation to the most general solution of the Renormalization Group Equation if 
gluinos and squarks are widely separated in mass. In section 3 we have presented the 
(formal) solution to this most general evolution, where each stage of the evolution 
produces shapes of the distributions which are used as input for the next regions. To 
be specific, let's consider the most general (Qo, Qi)ap - (Qi, Q2)sap - {Q2,Qf)ESAP 
supersymmetric evolution. We are allowed to vary the two matching scales Qi and 
Q2. Partons become massless -this is the content of a step- approximation to the gen- 
eral solution- as soon as we step into a new region. In this respect, we should remark 
that threshold effects, in order to be kept fully into account, would require some 
knowledge on the way supersymmetry is broken (or restored) as we raise the energy 
scale. This goes beyond our actual understanding of the theory. From a perturbative 
viewpoint, additional anomalous dimensions -needed to perform a matching between 
the various regions - may be needed, especially in moving from the SAP region to the 
ESAP region. Similarly to QCD, it is expected that in leading order these effects are 
negligible. In the analysis that follows we will consider the following scenarios: a) a 
light gluino and a heavy squark; b) a heavier gluino and a much heavier squark c) a 
heavy gluino and a heavier squark at an extremely large evolution scale. Being the 
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violations to scaling logarithmic, it is expected that the enhancements of the super- 
symmetric distributions will appear on an extremely large final evolution scale. As we 
are going to see, the enhancements -especially for gluinos- are larger at small-x and 
for smaller gluino masses. By varying the mass of the scalar quarks -here assumed to 
be degenerate, just for simplicity- we find that the modifications of the distributions 
are significant. We have kept nj = 6 in all the runs. Larger squark masses reduce 
drastically the maxima of the distributions, as illustrated below. 

In Fig. |15|we show the dependence of the gluino density using a general AP-SAP- 
ESAP evolution. We have 2 matching parameters: m2\ and rriq. Here we have varied 
both matching scales and the changes induced on the distributions are shown to be 
small. A more pronounced modification on the amplitudes can be observed from 
Fig. |T^, where we plot the non-singlet squark density for a light gluino and a heavier 
squark. The same pattern is observed in Fig. |1^ where we vary the mass of the squark 
at fixed gluino mass. The same effect of reduction of the distributions takes place as 
we raise the squark mass. In Fig. |1^ we compare 3 models of evolution: the general 
susy evolution, the partial susy evolution (AP-SAP) and the combined AP-ESAP 
evolution where we switch on an N=l evolution right on top of the regular AP run. 
It is observed that a run with partial susy generates distributions of gluinos which are 
smaller than the corresponding fully supersymmetric ones. It is also evident that a 
mixed AP-SAP-ESAP evolution and an AP-ESAP run generates similar results given 
the parameters chosen in this case (light gluino, heavier squarks). In Fig. ^ we plot 
the gluino distribution for a varying final evolution scale. The dependence is shown 
not to be significant. Fig. ^ and ^ show a similar pattern when we vary both the 
gluino mass and the squark mass in the general susy evolution. In these two figures 
we have chosen to vary both matching scales (gluino and squark masses), without 
any appreciable modification on the result. We conclude that the gluino distribution 
is not much affected by changes in the paramters of the evolution. The same type 
of numerical study (for an AP-SAP-ESAP run) is carried out in Fig. but for the 
singlet squark distribution. As we raise the squark mass, the distribution is lowered 
considerably. Here we have chosen a light gluino in the first range of the evolution. 
It is shown that the squark density is peaked at small-x and larger then a gluino 
density (light gluino) for the same x-value, at least as far as rusquark < 100 GeV. 



9 Conclusions 

Scaling violations induced by a supersymmetric evolution of the parton distributions 
have been studied in the context of a general supersymmetric evolution. The model 
distributions illustrated in this work have been generated within a radiative model, 
using the mixing of the QCD anomalous dimensions with the new states included in 
the = 1 as we cross each supersymmetric region in the evolution. Gluinos, and 
singlet squark distributions are strongly enhanced at small-x, have a faster decay at 



13 



larger-x compared to the gluon distribution and show dependence on the parameters 
(the two matching scales) of the evolution. The dependence is less pronounced for 
gluinos compared to squarks. The regular QCD distributions also show a small 
dependence on the parameters of the evolution, although these effects might be of 
minor phenomenological relevence at current energies. We will presente elsewhere 
applications of these results to the computation of supersymmetric cross sections. 
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11 Appendix 

The ESAP kernels are given by 



P(0) 
99 

p(0) _ 
p(0) _ 

p(0) _ 

p(0) _ 
^Xs — 

p(0)S ^ 
99 

p(0) _ 

^qX - 

p{0)S 
99 



\ +--2 + x(l-x) 



nx [2x - 1] 

C ' 



-1-X+ 



1-x) 
Ca(1 - x) 
Cf 

'x" + (1 - xf 



n/(l-x) 





= Cf[^ 


P(0) 
9S 


= Cf 


piO) 

sg 


= ^/[l 


p(0) 
^sX 


= UfX 


p(0) 
sg 


= Cfx 


P(0) 

ss 


= Cf( 



[1-x) 
2 



-S{l-x) 



+ 



(l-xU-'-'+-2^^'-''\ 



1 - (a;2 + (1 - 



(l-x). 



2 +CfS{1-x). 



(34) 



14 



The SAP kernels are given by 
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Figure 2: xG{x) for a regular AP versus an ESAP evolution, with m2\ = 20 GeV. 
The final evolution scale is Q/ = 500 GeV 




Figure 3: xuy{x) for a regular AP versus an AP-ESAP evolution, with m2x — 20 
GeV. The final evolution scale is Q/ = 500 GeV 
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Figure 4: xq~^{x) and xq^{x) for a regular AP versus an AP-ESAP evolution, with 
= 20 GeV. The final evolution scale is Qf — 500 GeV 
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Figure 5: xX{x) for the AP-ESAP evolution, with rriq — 20, 100, 200 GeV respectively. 
The final evolution scale is Q/ = 500 GeV 
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Figure 6: Non singlet distributions for the AP-ESAP evolution, with vriq = 20, 200 
GeV for xuv{x) and for the non-singlet distribution of scalar quarks {rriq — 200 GeV). 
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Figure 7: x in the AP-ESAP evolution, with = 20, 100, 200 GeV. The final 

evolution scale is Qf = 500 GeV 
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Figure 8: x q+{x) in the AP-ESAP evolution, with = 20, 100, 200 GeV. The final 
evolution scale is Qf = 500 GeV 
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Figure 9: xG{x) in the AP-ESAP evolution, with a squark mass rUg = 20 GeV. The 
final evolution scales are 1,2, and 5 TeV respectively 




Figure 10: xG{x) in the AP-ESAP evolution, with a squark mass rriq = 200 GeV. 
The final evolution scales are 1,2, and 5 TeV respectively 
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Figure 11: xX{x) in the AP-ESAP evolution, with a squark mass rrig = 200 GeV. 
The final evolution scales are 1,2, and 5 TeV respectively 
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Figure 12: xq (x) in the AP-ESAP evolution, with a squark mass = 20 GeV. 
The final evolution scales are 1,2, and 5 TeV respectively 
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Figure 13: xq~{x) in the AP-ESAP evolution, with a squark mass rriq = 20 GeV. 
The final evolution scales are 1,2, and 5 TeV respectively 




Figure 14: xf{x) in the AP-ESAP evolution with a very large final scale Qf = 10^ 
GeV and with a squark mass rriq = 100 GeV. Shown are the non-singlet squark, the 
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Figure 15: xX{x) using a general evolution. The matching parameters are denoted 
by two entries {m2x,mg) (in GeV) and Qj = 500 GeV 
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Figure 16: xq (x) using a general evolution. The matching parameters are indicated 
as in Fig. 1 and Qf = 500 GeV 
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Figure 17: Dependence of the singlet squark distribution on the squark mass for 
msA = 20 GeV and rrig = 70,90,200,300 and 400 GeV. We have chosen Q/ = 500 
TeV 
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Figure 18: 3 gluino distributions according to the 3 possible evolution models, shown 
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Figure 19: Gluino distribution xX{x) for = 40 GeV and nig = 100 GeV and 
Qj = 500, 800 Gev and 2 TeV respectively 




Figure 20: 3 gluino distributions according to the 3 possible evolution models, shown 
for m2x = 100 and 20 GeV and for m5=300,400 GeV. We have chosen Qf = 1 TeV 




Figure 21: 3 gluino distributions according to the 3 possible evolution models, shown 
for m2x = 20 GeV and rUg = 70, 90 and 200 GeV. We have chosen Q/ = 500 TeV 
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Figure 22: x q~^{x) distributions for m2\ = 40 GeV and = 70,90,200 and 
400 GeV respectively. Shown is also a gluino xX{x) distribution. We have chosen 
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